AD-A150  033 


UNIVERSITY  OF  NEVADA  -  RENO  -  NEVADA  -  89557 


‘VjtO' 


Mackoy  School  of  Mines 
Sei smalogi cal  Lahorotary 


ANNUAL 

TECHNICAL  REPORT  NO.  1 


Telephone  (702)  784-4975 


03  December  1979 


ARPA  Order  No . 


Program  Code 


Name  of  Contractor 


University  of  Nevada,  Reno 


Effective  Date  of  Contract 


01  Oct  1978 


Contract  Expiration  Date 


30  Sep  1979 


Amount  of  Contract  Dollars 


$39,226 


Contract  Number 


F49620-77-C-0070 


Principal  Investigator 


William  A.  Peppin,  702-734-4975 


Program  Manager 


William  A.  Peppin 


Short  Title  of  Work 


Close-in  and  Regional  Source 
Studies  for  Seismic  Discrimination 


Sponsored  by 


Advanced  Research  Projects  Agency  (D0D) 


ARPA  Order  No.  3291 


Monitored  by  AF0SR  under  Contract  //  F49620-77-C-0070 


*4J 

ur 


AP^-P  ' 


i  ■,  i).  n'lwVii  i  iiD  (A) 


s 


DTIC 

ELECTE 

JAN  2  4  1985  I 


85  oi  15  049 


A  DIVISION  OF  THE  UNIVERSITY  OF  NEVADA  SYSTEM 


"  L  >  V  »  -  r-  ann  r  i*Trr* 


TABLE  OF  CONTENTS 


Page 


TECHNICAL  REPORT  SUMMARY 

9 

1.  CO/E  FOR  SOURCE  INVERSION 

2.  DEVELOPMENT  ON  THE  SEISMOLOGICAL  LABORATORY 

COMPUTER  SYSTEM 

3.  SEISMIC  DISCRIMINATION  AT  NEAR-REGIONAL  DISTANCES 

3.1  PN  versus  Rayleigh  Discrimination  at  Jamestown 
3*2  P  versus  S  Discriminant  at  Jamestown 

3.3  Seismic  Discrimination  at  Jamestown:  Summary 

4.  SEISMIC  DISCRIMINATION  AT  REGIONAL  DISTANCES  USING  PN 

5.  ANALYSIS  OF  CLOSE-IN  ACCELEROMETER  RECORDS  OF 

MEGATON  EXPLOSIONS 

5.1  Source  in  a  Layer 

5.2  Source  Below  a  Single  Layer 

5.3  Helmberger-Hadley  Model 

5.4  Extended  Helmberger-Hadley  Model 

6.  PN  DISCRIMINATION  IN  A  SHIELDLIKE  ENVIRONMENT 

7.  PERIPHERAL  STUDIES  OF  INTEREST  TO  THIS  CONTRACT 

7.1  Digital  Data  Acquisition  and  Spectral  Corner  Frequencies 

7.2  Analysis  of  Explosion  Surface  Waves  for  Phase  Velocity 
REFERENCES  CITED 


f  Co»r  ' 

V  • 


Accession  For 

NT  IS  CFA&I 
DTIC  TAB 
Unpnnounced 
Justif ication_ 


XT 

□ 

a 


By- 


Distribution/ 
Availability  Codes 
Avail  and/or 
Special 


1 

2 

4 

5 
5 
8 

11 

13 

16 

17 

18 
19 
19 
21 
22 
22 

24 

25 


BEST 

AVAILABLE  COPY 


TECHNICAL  REPORT  SUMMARY 


This  contract  has  as  Its  primary  objectives  the  Investigation  of  close-in 
and  near-regional  <1.e.,  0  to  600  km>  seismic  data  for  the  purpose  of  Implement¬ 
ing  and  understanding  seismic  discriminants .  by  which  we  mean  any  methods  based 
on  the  analysis  of  seismic  waves  which  permit  unambiguous  Identification  of  se¬ 
ismic  explosions  as  either  underground  explosions  or  earthquakes.  Efforts  dur¬ 
ing  this  contract  period  have  emphasized  the  development  and  testing  of  a  n' 
and  rather  sophisticated  Inversion  method  for  seismic  sources,  the  seismic  mo¬ 
ment  tensor  method,  whose  use  was  first  attempted  by  Stump  and  Johnson  (1977)  at 
these  distance  ranges.  Material  presented  Includes  the  following.  In  Section  1 
we  present  comments  about  the  development  of  a  major  new  computer  program  spec¬ 
ifically  built  to  produce  the  more  realistic  Green's  functions  needed  for  moment 
tensor  analysis.  Section  2  Is  a  discussion  of  developments  accomplished  on  the 
Selsmologlcal  Laboratory  computer  system  In  support  of  this  contract.  In  Sec¬ 
tion  3  Is  presented  a  discussion  of  two  near-regional  discriminants  based  on  new 
wideband  data*  the  PN-Raylelgh  discriminant  of  McEvIlly  and  Peppln  (1972),  and 
an  S  to  P  ratio  method  that  has  been  Investigated  under  different  names  by  sev¬ 
eral  authors.  In  Section  4  Is  described  an  attempt  to  effect  seismic  discrimi¬ 
nation  through  use  only  of  the  PN  phase  using  the  wideband  three-component  data 
of  the  Lawrence  Livermore  Laboratory.  In  Section  6  Green's  functions  for  moment 
tensor  analysis  of  the  close-in  JORUM  and  HANDLEY  accelerometer  data  are  dis¬ 
cussed.  This  work  extends  the  previous  work  of  the  author  and  a  number  of  oth¬ 
ers  on  this  Important  data  set.  In  Section  6  we  discuss  efforts  directed  at 
testing  the  PN  discriminant  described  In  Section  4  over  a  "shield! Ike*  structure 
based  on  large  underground  explosions  In  Missouri.  In  Section  7  we  discuss  some 
work  that  Is  peripheral,  but  related  to  this  contract  Involving  the  study  of  the 
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earthquake  source  mechanism,  Including  the  recent  series  of  earthquakes  along 
the  eastern  Sierra  front  and  Rayleigh  wave  dispersion  on  the  test  site. 

1.  CODE  FOR  SOURCE  INVERSION 

Midway  through  this  fiscal  year  I  made  operational  a  code  to  compute  ground 
displacement  from  general  seismic  sources.  This  code,  named  MEXEC,  Is  a  full 
general Izatlon  of  the  earlier  one  used  for  the  computation  of  ground  motion  from 
explosion  sources.  The  output  Is  directly  compatible  with  Brian  Stump's  moment 
tensor  Inversion  code  (Stump  and  Johnson,  1977}  Stump,  1979).  The  general  phi¬ 
losophy  of  this  approach  to  seismic  sources  Is  distinct  from  the  more  commonly- 
applied  ones.  Usually  the  assumption  Is  made  that  the  source  Is  either  a  pure 
explosion  or  a  pure  earthquake  with  fault  planes  having  some  orientation;  then 
some  method  such  as  generalized  ray  theory  Is  used  to  propagate  energy  from  this 
source  to  the  receiver  (see  the  treatment  by  Johnson  and  McEvIlly,  1973  for  ex¬ 
ample).  In  moment  tensor  analysis,  we  make  no  such  assumption.  Rather,  appeal¬ 
ing  to  the  theory  of  Burrldge  and  Knopoff  (1964),  we  assume  the  source  to  be  re¬ 
presentable  by  a  system  of  burled  body  forces.  These  body  forces  are  expanded 
about  a  convenient  point  In  a  series  known  as  the  moment  tensor  representat Ion 
(Stump  and  Johnson,  1977).  Then,  observed  seismograms  are  fitted  by  such  a  li¬ 
near  combination  of  these  force  moments  that  the  residual  of  predicted  from  ob¬ 
served  seismograms  Is  minimized  In  a  least-squares  sense.  This  procedure  ef¬ 
fects  an  Inversion  for  the  components  of  the  “seismic  moment  tensor",  which  for 
small  sources  can  be  taken  as  particular  linear  combinations  of  the  derivatives 
of  the  Green's  function  for  the  problem.  What  Is  significant  about  this  proce¬ 
dure  Is  that  no  a  priori  assumption  about  the  source  Is  made  with  regard  to  Its 
Identity  as  an  explosion  or  earthquake!  that  Information  emerges  from  the 


PAGE  3 


Inversion  together  with  error  estimates  as  to  how  precisely  the  moment  compo¬ 
nents  are  known.  In  other  words*  Inversion  Is  done  and  we  can  say  that  a  cer¬ 
tain  percent  of  the  source  (with  some  uncertainty)  Is  an  explosion,  some  other 
certain  percent  Is  contained  In,  e.g.,  a  dip-slip  earthquake,  and  so  forth. 
This  Is  a  statement  of  the  central  problem  of  seismic  dl scrlmlnatlon,  and  so  Is 
the  natural  method  to  use.  To  date,  no  work  on  seismic  discrimination  has  been 
done  at  the  close-in  and  near-regional  distance  ranges  using  the  moment  tenser 
method.  The  MEXEC  code  and  a  host  of  data  collected  by  this  author  have  been 
brought  together  In  an  effort  to  test  thoroughly  the  efficiency  of  the  method 
and,  It  Is  hoped,  gain  some  new  Insight  Into  the  nature  of  the  seismic  source. 

In  the  course  of  writing  the  MEXEC  code,  I  have  redesigned  the  parts  where 
the  source  terms  are  entered  In  such  a  way  that  computation  of  all  of  the 
higher-order  moment  source  types  Is  trivial.  Although  only  those  of  first  order 
are  computed  In  the  present  version  of  the  code.  It  Is  a  matter  of  changing  a 
few  lines  of  FORTRAN  In  a  single  subroutine  to  treat  any  of  these  components.  A 
considerable  time  during  this  contract  period  was  also  spent  In  documenting  the 
code  and  some  of  the  more  sophisticated  numerical  difficulties  that  have  come 
up.  As  a  result  of  this  effort,  I  believe  that  the  code  Is  finally  working  cor¬ 
rectly.  Some  quite  critical  tests  have  been  made  although  I  cannot  yet  certify 
the  correctness  of  the  code.  A  central  problem  In  bringing  up  this  code  was  a 
collection  of  quite  difficult  numerical  problems  that  had  to  be  overcome.  The 
Cagnlard-deHoop  method  permits  an  elegant  and  simple  way  to  write  down  the  solu¬ 
tions  to  difficult  problems  In  wave  propagation;  however,  transforming  this  de¬ 
ceptively  simple  expressions  Into  an  operation  computer  code  Is  altogether 
another  matter  (the  solutions  were  developed  In  about  300  hours  of  work,  but  the 
program  has  consumed  no  less  than  5,000  hours  of  development  time). 
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I  should  emphasize  why  I  endured  such  a  long  development  effort  when  gener¬ 
alized-ray  codes  were  available  from  groups  at  Caltech  and  elsewhere.  The 
answer  Is  twofold*  first,  I  never  found  the  documentation  provided  by  the  Cal¬ 
tech  croup  sufficient  to  understand  completely  what  was  being  done,  so  I  would 
have  had  great  difficulty  understanding  their  codes;  second,  my  code  differs 
from  the  more  usual  ones  Involving  genera  1 1  zed-ray  analysis  In  that  It  Is  exact 
at  all  distance  ranges,  no  approximation  having  been  made  anywhere  In  the  ana¬ 
lysis.  This  permits  Investigations  of  seismic  records  at  very  close  distances 
to  the  source  without  the  worry  of  having  to  know  which  terms  In  the  asymptotic 
series  to  retain. 

Cooperative  studies  are  underway  with  Lawrence  Livermore  Laboratory,  where 
I  will  again  be  a  consultant  for  the  time  period  of  this  contract.  MEXEC  has 
bean  placed  on  LTSS  at  Livermore,  and  on  the  Prime  computer  In  the  Geophysics 
division  there.  It  has  also  been  made  operational  at  Klrtland  Air  Force  Base 
for  use  by  Brian  Stump  for  projects  Involved  In  his  reponslbl  1  It les  and  for  mo¬ 
ment  tensor  analysis.  In  Sections  4  and  6  are  presented  examples  of  output  from 
the  MEXEC  code. 

2.  DEVELOPMENTS  ON  THE  SEISMOLOGICAL  LABORATORY  COMPUTER  SYSTEM 

During  this  contract  period  a  considerable  body  of  FORTRAN  software  was 
created  for  digital  signal  processing  on  the  Selsmologlcal  Laboratory  POP  11/34 
computer.  This  software  was  developed  In  the  course  of  the  completion  of  three 
Master's  degrees  by  grad  students,  and  for  the  processing  of  the  data  from  the 
digital  event  recorders  developed  with  support  from  AFOSR.  Under  Keith  Priest¬ 
ley's  contract  with  AFOSR,  codes  were  developed  for  determining  the  structure  of 
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the  Basin  and  Range  using  surface-wave  analysis,  gravity,  and  seismic  refraction 
profiles.  In  support  of  this  contract,  I  have  developed  the  software  for  pro¬ 
cessing  records  obtained  from  the  Berkeley  tape  system  (see  Section  4).  The 
system  has  been  upgraded  by  acquisition  of  a  maintenance  contract  from  Digital 
Equipment,  by  Implementation  of  the  new  release  of  RSX11-M,  purchase  of  a  fast 
compiler  FORTRAN  IV  PLUS,  and  a  remote  hardcopy  device  for  the  Tektronix  graph¬ 
ics  terminal.  The  computer  has  served  all  of  our  real-time  and  other  processing 
needs  very  well,  and  we  anticipate  no  significant  need  of  additional  hardware 
and  software  support  for  the  objectives  of  this  contract.  During  the  contract 
period  our  computer  systems  analyst  position  changed  hands  from  Ron  Sheen  to 
Dennis  Ghlgllerl,  who  has  proven  an  able  replacement.  The  staff  uses  the  com¬ 
puter  almost  dally,  a  clear  testimonial  to  the  simplicity  and  usefulness  of  the 
machine. 


3.  SEISMIC  DISCRIMINATION  AT  NEAR-REGIONAL  DISTANCES 


3.1.  PN  versus  Rayleigh  Discrimination  at  Jamestown 


McEvIlly  and  Peppln  (1972)  and  Peppln  and  Me  Evilly  (1973)  found  a  success¬ 
ful  seismic  discriminant  between  underground  nuclear  explosions  at  Nevada  Test 
Site  (NTS)  and  earthquakes.  Making  use  of  the  wideband  systems  at  Berkeley, 
Mina ,  and  Elko,  Nevada,  Kanab,  Utah,  and  Landers,  California  (distance  range  200 
to  660  km  from  NTS),  almost  all  events  were  successfully  separated  by  a  good 
margin,  providing  one  of  the  few  effective  near-regional  discriminants.  The 
discriminant  compared  the  maximum  amplitude  of  the  PN  phase  In  a  pass  band  0.6 
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to  5  Hz  with  maximum  filtered  (0.02  to  0.01  Hz)  Rayleigh  wave  motion.  Surface 
waves  were  detectable  for  explosions  down  to  a  body-wave  magnitude  of  about  3.6, 
and  for  earthquakes  somewhat  lower. 

Now  the  discriminant  was  found  to  be  effective  at  Berkeley  and  Mina,  both 
northwest  of  the  test  site,  but  less  effective  In  the  azimuths  towards  Kanab  and 
Landers  (l.e.,  due  east  and  due  south  from  the  test  site).  It  Is  of  Interest  to 
know  If  this  discriminant  works  In  general.  Therefore,  when  the  University  of 
California  at  Berkeley  commenced  operation  of  a  wideband  system  at  Jamestown, 
also  NW  of  the  test  site  and  In  the  Sierra  foothills,  an  opportunity  was  provid¬ 
ed  to  study  the  discriminant  at  yet  another  wideband  station.  By  this  time  suf¬ 
ficient  data  has  been  recorded  on  the  Jamestown  system  to  test  the  PN-Raylelgh 
discriminant  rather  thoroughly.  It  Is  noteworthy  that  the  Jamestown  station 
simulated  an  SRO  site,  because  the  data,  although  recorded  analog,  are  kept  for 
a  variety  of  gains  and  bandwldths.  Thus,  data  Is  available  for  earthquakes  any¬ 
where  In  the  western  United  States  In  the  magnitude  ranges  from  about  3  to  6.6. 
The  pass  bands  used  for  the  analysis  described  here  were  the  hlghgaln  shortperl- 

H 

od  (0.2  to  1.0  /k,  "SPZ" ),  wideband  velocity  (0.026  to  10  Hz,  0.6 
volts/mlcron/sec,  "BBV“),  and  wideband  displacement  taken  at  two  gain  levels 
(0.026  to  10  Hz,  0.6  or  0.0026  volts/micron,  "DHG"  or  "DLG"  for  high  and  low 
gains,  respectl vely).  Whole-record  data  for  explosions  In  the  local  magnitude 
range  3.0  to  6.6  Is  available  for  the  continuously-recorded  data  set  running 
frcm  1676  to  the  present  time. 

From  Jun  1975  to  Feb  1979  IJ  events  were  selected  for  analysis.  As  In  the 
previous  studies,  the  primary  difficulty  Is  In  finding  suitably-recorded  earth¬ 
quakes  on  the  test  site.  There  was  none  at  all  on  the  test  site  except  hole 


PAGE  7 


collapses.  Consequently!  following  the  method  of  McEvIlly  and  Peppln  (1972),  a 
search  was  made  for  events  at  around  400  km  from  Jamestown.  Earthquakes  found 
spanned  a  large  range  of  azimuths  from  NW  (near  Shasta  Dam,  California)  through 
northern  and  central  Nevada,  Death  Valley,  China  Lake,  and  the  San  Fernando  Val¬ 
ley  (about  225  degree  range  of  azimuths).  Of  the  events  selected  48  had  PN  and 
Rayleigh  wave  signals  that  could  be  studied  (Table  1).  The  smallest  magnitudes 
were  slightly  below  ML  4.0  (which  Is  not  as  low  as  deemed  desirable  for  seismic 
verification  purposes). 

Figure  3.1  Is  a  plot  similar  to  those  presented  by  McEvIlly  and  Peppln 
(1972)  and  Peppin  and  McEvIlly  (1973).  Numerals  by  each  symbol  are  event 
numbers  In  Table  1.  Surface-wave  amplitudes  were  read  on  either  the  BBV,  DHG, 
or  DLG  channels  off  vlslcorder  playouts  fom  analog  tape  after  filtering  (0.02  to 
.10  Hz,  48  db/octave  rolloff).  Presented  for  a  Rayleigh  wave  measurement  Is  the 
maximum  peak-to-peak  amplitude  on  the  filtered  trace  of  the  DHG  channel.  For 
s;re  cases  (asterisks  In  Table  1)  It  was  Impossible  to  obtain  a  reading  on  the 
DHG  channel.  In  that  case  the  reading  was  obtained  either  from  DLG  through  the 
static  gain  adjustment  between  them  of  200,  or  through  an  empirical  adjustment 

DHG  -  BBV/( 11 . 1  ±  2.6) 

obtained  for  the  22  events  having  a  useable  surface  wave  on  both  the  DHG  and  BBV 
channels.  Presented  for  the  PN  measurement  Is  the  maximum  peak-to-peak  ampli¬ 
tude  determined  on  SP2,  filtered  0.5  to  3.0  Hz  (24  db/octave  rolloff).  The 
larger  events  clipped  on  this  trace,  so  an  empirical  adjustment  was  made  from 
the  BBV  and  DHG  traces  as 


LOG  RRTLEIGH  RMP 


0  12  3  4  5 

LOG  PN  AMPLITUDE 

Figure  3.1 .  Pn  versus  Rayleigh  discrimination  for  events  about  400  km 
from  the  Berkeley  station  Jamestown.  Event  numbers  refer  to  Table  1. 


IDENTIFIER  EVENT  NO 

PN  AMP 

RAYLEIGH  AMP 

COMMENTS 

1975  JUN  27  0727 

1 

32 

88* 

COLLAPSE 

1975  JUN  28  1948 

2 

28 

77* 

COLLAPSE 

1975  JUL  01  0451 

3 

18 

30 

COLLAPSE 

1975  JUL  01  1814 

4 

40 

248* 

COLLAPSE 

1976  SEP  06  1700 

5 

104 

18 

EXPLOSION 

1975  OCT  24  1712 

6 

323* 

57 

EXPLOSION 

1975  OCT  28  1430 

7 

17000* 

16800 

EXPLOSION 

1975  NOV  20  1600 

8 

2450 

3440 

EXPLOSION 

1975  NOV  26  1530 

9 

104 

47 

EXPLOSION 

1976  JAN  03  1915 

10 

6420* 

25600 

EXPLOSION 

1976  JAN  04  0118 

11 

42 

76 

COLLAPSE 

1976  JAN  04  1616 

12 

18+ 

60+ 

COLLAPSE 

1976  JAN  17  2140 

13 

22 

40 

COLLAPSE 

1976  FEB  04  1420 

14 

1275* 

1160* 

EXPLOSION 

1976  FE3  04  1440 

15 

1500 

1309* 

EXPLOSION 

1976  FEB  12  1445 

16 

6900 

5525 

EXPLOSION 

1976  FEB  14  1130 

17 

3434* 

24000 

EXPLOSION 

1976  FEB  19  1701 

18 

28 

10 

EARTHQUAKE, 

A2 

130,  MAG  4.3,  D  332 

1976  MAR  09  1400 

19 

3485* 

6400 

EXPLOSION 

1976  MAR  14  1430 

20 

7480* 

17600 

EXPLOSION 

1976  MAR  14  1525 

21 

56 

200 

COLLAPSE 

1976  MAY  12  1950 

22 

144 

88 

EXPLOSION 

1976  JUN  07  0035 

23 

2 

44 

EARTHQUAKE, 

AZ 

115,  MAG  4.1,  D  392 

1976  AUG  26  1430 

24 

408 

248 

EXPLOSION 

1976  DEC  08  1450 

25 

136 

112 

EXPLOSION 

1976  OEC  21  1509 

26 

104 

1 

EXPLOSION 

1977  JAN  13  0719 

27 

74 

10* 

EARTHQUAKE, 

AZ 

170,  MAG  3.8,  D  370 

1977  JAN  13  2009 

28 

44+ 

28 

EARTHQUAKE, 

AZ 

330,  MAG  3.7,  D  377 

1977  FEB  16  1753 

29 

96 

22* 

EXPLOSION 

1977  AUG  04  1640 

30 

328 

440 

EXPLOSION 

1977  AUG  12  0220 

31 

224 

180 

EARTHQUAKE , 

AZ 

170,  MAG  4.8,  D  421 

1977  AUG  19  1755 

32 

744 

1386* 

EXPLOSION 

1977  SEP  15  1436 

33 

72 

11* 

EXPLOSION 

1977  OCT  26  1415 

34 

56 

76 

EXPLOSION 

1977  NOV  10  0236 

35 

48+ 

15 

EARTHQUAKE, 

AZ 

126,  MAG  4.0,  D  302 

1977  NOV  17  1930 

36 

88 

26 

EXPLOSION 

1977  DEC  14  1530 

37 

3264* 

1400 

EXPLOSION 

1978  FEB  14  0435 

38 

328 

256 

EARTHQUAKE , 

AZ 

45,  MAG  4.8,  D  340 

1978  MAY  23  0547 

39 

104 

74 

EARTHQUAKE , 

AZ 

20,  MAG  4.6,  D  434 

1978  JUN  16  0421 

40 

64 

15 

EARTHQUAKE , 

AZ 

170,  MAG  4.3,  D  346 

1978  AUG  01  0902 

41 

56 

700 

EARTHQUAKE , 

AZ 

350,  MAG  4.6,  D  397 

1978  AUG  01  0947 

42 

30 

528 

EARTHQUAKE, 

AZ  350,  MAG  4.5,  D  397 

1978  AUG  01  1026 

43 

44 

58 

EARTHQUAKE , 

AZ 

350,  MAG  4.2,  D  397 

1978  AUG  31  1400 

44 

5440* 

1400 

EXPLOSION 

1978  AUG  31  2356 

45 

42 

86 

COLLAPSE 

1978  NOV  02  1625 

46 

168 

55 

EXPLOSION 

1979  JAN  06  0120 

47 

80+ 

39 

EARTHQUAKE , 

AZ 

75,  MAG  4.2,  D  377 

1979  JAN  24  1800 

48 

200 

66  * 

EXPLOSION 

*  -  READINGS  BY  CONVERSION  TO  STANDARP  ONES  (SEE  TEXT) 
+  -  PCCR  READINGS 


Table  1.  Basic  data  taken  for  the  investigation  of  the  Pn  versus 
Rayleigh  discriminant.  "AZ\  ;,MAG" ,  AND  "D"  are  azimuth  from  James¬ 
town  to  the  event,  the  local  magnitude,  and  distance  in  km. 
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SPZ  -  (17.1  t  6)BBV  (IB  cases), 

or 

SPZ  -  (34.3  ±  12 )DHG  (8  cases). 

The  observations  starred  In  Table  1  are  those  for  which  adjusted  SPZ  read¬ 
ings  are  given  (smaller  symbols  In  the  figure). 

The  results  confirm  my  suspicion  of  many  years'  standing:  the  discriminant 
falls.  As  a  group  only  hole  collapses  are  separated  from  the  explosions.  This 
highlights  the  discussion  presented  In  Peppln  (1974).  We  are  faced  with  an  em¬ 
pirical  method  that  seems  to  Imply  a  possible  discrimination  of  earthquakes  from 
explosions;  but  we  are  not  sure  sure  that  the  discriminant  will  be  effect. ve  If 
applied  to  a  region  of  Interest  elsewhere  (l.e.,  the  Soviet  Union).  Only  this 
year  hive  I  obtained  a  means  to  Investigate  the  theoretical  causes  for  the 
PN-Raylelgh  discriminant  (see  next  section).  However,  given  this  negative  out¬ 
come,  the  results  of  that  analysis  would  appear  of  less  Interest. 

Note  that  the  situation  we  have  presented  Is  really  bad  for  the  effective¬ 
ness  of  this  discriminant,  because  all  of  the  stations  Berkeley,  Jamestown,  and 
Mina  are  roughly  the  same  azimuth  from  the  test  site.  The  fact  that  the  near 
and  far  stations  seem  to  provide  effective  discrimination  and  the  middle  one  not 
Is  very  bad  news  for  him  desirous  of  applying  this  discriminant  to  any  case  of 
real  Interest. 

3.2.  P  versus  S  Discriminant  at  Jamestown 
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In  performing  the  data  search  for  the  previous  section,  I  noticed  that  one 
could  tell  at  a  glance  which  events  were  earthquakes  and  which  were  not.  In 
virtually  all  cases  the  earthquakes  had  visibly  and  obviously  larger  amplitudes 
of  the  Sg  waves  -Ive  tu  Pg.  Various  authors  have  commented  on  this  begin¬ 
ning  with  the  first  paper  on  discrimination  (Leet,  1962)*  since  explosions  are 
symmetric  sources,  they  should  generate  less  S  energy  than  a  shear  source. 
Indeed,  this  Is  taken  to  be  the  cause  of  the  MsSmb  discriminant,  since  both  P 
and  S  leaving  the  source  generate  Rayleigh  waves  (Douglas  et.  al.,  l'37t). 
Murphy  and  Lahoud  (1976)  have  given  a  more  detailed  account  of  this  phenomenon 
for  near-regional  records  of  NTS  earthquakes  and  explosions. 

Motivated  by  the  above,  I  made  a  comparison  of  the  ratio  of  Sg  to  Pg  maxi¬ 
mum  amplitude  for  a  superset  of  the  events  used  In  Section  2.1.  The  discrimina¬ 
tion  technique  I  describe  below  Is  simple  to  apply  among  Its  advantages,  and  I 
have  been  able  to  place  quantitative  bounds  on  Its  range  of  effectiveness. 

In  a  study  of  60  events  1  first  attempted  to  see  If  a  spectral  estimate 
would  be  better  than  an  analog  one.  However,  visually  It  appeared  that  discrim¬ 
ination  based  on  wideband  spectral  power  would  be  less  effective  than  a  measure¬ 
ment  of  maximum  amplitude  (and  Indeed  Peppln  and  McEvIlly,  1973  had  found  no 
luck  In  their  attempts  to  develop  a  discriminant  on  wideband  spectral  averages 
of  the  Pg  phase). 

The  scheme  adopted  was  as  follows.  First,  draw  an  envelope  around  the 
bursts  of  Pg  and  Sg  energy  (Figure  3.2);  this  Is  an  attempt  to  smooth  the  data 
partially  rather  than  measure  the  Isolated,  sharp  maximum  peaks.  Measure  the 
0-pk  height  of  the  envelope  at  the  onset  times  of  the  Pg  and  Sg  phases  (l.e., 


P»\“-  w*»  • 


-v- 


the  ratio  B  to  A  In  the  Figure)  as  seen  on  the  SPZ  trace  played  out  on  a  vlsl- 
corder.  Plot  this  ratio  as  a  function  of  local  magnitude.  We  are  dealing  with 
a  tremendous  range  of  amplitudes,  so  empirical  adjustments  are  made  to  connect 
the  different  components  recorded  by  the  Jamestown  system  as  In  the  previous 
section.  These  adjustments  were 

SPZ  ■  DHG(17.6  ±  7.8)  for  P,  23  observations 
SPZ  ■  DHG(13,6  ±  6.1)  for  S,  23  observations 
SPZ  M  BBV<13.9  ±  4.2)  for  P,  14  observations 
SPZ  *  BBV(13.0  ±  4.5)  for  S,  24  observations 

Results  of  the  analysis  are  given  In  Table  2  and  Figure  3,\  Note  that  the 
explosions  separate  well  from  the  natural  earthquakes,  but  that  the  hole  col¬ 
lapses  extend  far  up  Into  the  earthquake  population.  All  explosions  show  an  S 
to  P  ratio  less  than  2  and  all  earthquakes  except  one  <  28  In  Table  2)  show  a 
ratio  as  least  as  great  as  2. 

In  order  to  assign  statistical  significance  to  the  degree  of  separation,  I 
computed  an  uncertainty  for  those  events  recorded  on  each  of  BBZ,  DHG,  or  SPZ. 
These  uncertainties  are  shown  In  Figure  3.4.  I  have  displaced  these  96*  error 
bars  left  or  right  of  the  event  In  question  for  visual  clarity  of  presentation. 
Standard  errors  amount  to  32*,  17%,  and  16*  of  the  mean  for  collapses,  explo¬ 
sions,  and  earthquakes,  respect  1 ve ly ,  where  we  have  Ignored  the  earthquake  with 
huge  error  bars  at  Shasta  that  Is  obviously  not  an  explosion  from  the  appearance 
of  the  records  taken.  In  Figures  3.3  and  3.4  we  have  drawn  a  horizontal  line 
separating  the  populations  at  a  ratio  of  2  and  96*  error  bars  20*  of  the  mean 
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around  either  side  of  It.  Of  sixty  events  ranging  In  magnitude  from  3.5  to  6.5 
13  are  definitely  Identified  as  not  explosions,  36  are  definitely  Identified  as 
explosions,  and  the  remaining  11  are  unresolved.  A  single  earthquake  <  28  In 
Table  2  and  Figure  3.4)  falls  Into  the  explosion  population.  However,  this 
e'/ent  (  23  In  Table  1  and  Figure  3.1)  Is  definitely  identified  as  an  earthquake 
by  the  PN-Raylelgh  discriminant.  For  discrimination  purposes,  it  is  not  a  fatal 
flaw  that  many  collapses  are  mlsldentlf led  as  explosions,  since  the  presence  of 
a  collapse  implies  (with  about  99%  probability)  an  explosion  within  the  last 
week.  Many  of  the  collapses  can  be  Identified  by  their  large  surface  wave  exci¬ 
tation  (Figure  3.1),  and  spectral  ratio  techniques  could  identify  many  more 
(most  collapses  show  obviously  longer  periods  in  their  Pg  waves  than  are  seen  In 
explosions  of  comparable  magnitude!  see  the  work  of  Stump,  1979  for  example). 


3.3.  Seismic  Discrimination  at  Jamestown!  Summary 

If  we  use  as  a  primary  discriminant  the  S/P  ratio  method  described  In  Sec¬ 
tion  3.2  and  The  PN-Raylelgh  discriminant  as  a  secondary  one,  28  of  30  explo¬ 
sions  are  Identified  with  certainty?  no  earthquake  in  a  region  extending  from 
northern  California,  through  Nevada,  and  back  into  southern  California  Is  Iden¬ 
tified  as  an  explosion.  Hole  collapses  remain  ambiguously  or  Incorrectly  iden¬ 
tified,  but  the  Identification  of  these  should  be  no  problem.  The  S/P  ratio 
method  has  some  compelling  and  obvious  advantages*  (1)  it  is  effective  to  mag¬ 
nitudes  as  small  as  can  be  seen  on  high-gain,  short-period  seismographs  at 
near-regional  distances,  l.e.,  local  magnitude  3.0  or  less;  (2)  It  Is  simple  to 
apply;  (3)  It  Is  reasonably  easy  to  understand  how  It  works;  (4)  plots  like 
Figure  3.3  can  be  carried  over  to  another  region  without  change  as  a  ratio  Is 


*  EXPLOSION 
O  COLLAPSE 
+  EARTHQUAKE 


COMMENTS 


IDENTIFIER  EVENT  NO  P  AMP  S  TO  P  RATIO 


1975 

JUN 

27 

0727 

1 

15. 

1.12* 

COLLAPSE 

1976 

JUN 

28 

0948 

2 

14. 

1.10* 

COLLAPSE 

1976 

JUL 

01 

0451 

3 

12. 

1.05 

COLLAPSE 

1976 

JUL 

01 

1814 

4 

39. 

2.68 

COLLAPSE 

1976 

SEP 

06 

1700 

5 

22. 

1.11 

EXPLOSION 

1976 

OCT 

24 

1712 

6 

67. 

1.24 

EXPLOSION 

1976 

OCT 

28 

1430 

7 

76. 

0.32 

EXPLOSION 

1976 

NOV 

18 

1530 

8 

12. 

1.75 

EXPLOSION 

1976 

NOV 

20 

1500 

9 

1400. 

0.50 

EXPLOSION 

1976 

NOV 

26 

1630 

10 

36. 

1.47 

EXPLOSION 

1976 

JAN 

03 

1915 

11 

4200. 

0.48 

EXPLOSION 

1976 

JAN 

04 

0118 

12 

13. 

1.26 

COLLAPSE 

1976 

JAN 

04 

1616 

13 

8. 

2.40 

COLLAPSE 

1976 

JAN 

18 

0720 

14 

14. 

1.75 

COLLAPSE 

1976 

JAN 

21 

1840 

15 

9. 

2.11 

COLLAPSE 

1976 

FEB 

03 

0015 

16 

8. 

1.33 

COLLAPSE 

1976 

FEB 

04 

1420 

17 

480. 

1.20 

EXPLOSION 

1976 

FEZ 

o4 

1440 

18 

416. 

1.58 

EXPLOSION 

1976 

FEB 

12 

1445 

19 

1400. 

0.71 

EXPLOSION 

1976 

FEB 

14 

1130 

20 

3800. 

0.53 

EXPLOSION 

1976 

FEB 

19 

1701 

21 

15. 

3.97 

EARTHQUAKE, 

AZ 

130,  MAG  4.3,  D  332 

1976 

FEB 

26 

1450 

22 

8. 

1.20 

EXPLOSION 

1976 

MAR 

09 

1400 

23 

2000. 

1.33 

EXPLOSION 

1976 

MAR 

09 

1655 

24 

16. 

1.26 

COLLAPSE 

1976 

MAR 

14 

1450 

25 

4400. 

0.45 

EXPLOSION 

1976 

MAR 

14 

1625 

26 

30. 

1.22* 

COLLAPSE 

1976 

MAY 

12 

1960 

27 

53. 

1.18 

EXPLOSION 

1976 

JUN 

07 

0036 

28 

13. 

1.77 

EARTHQUAKE , 

AZ 

115,  MAG  4.1,  D  392 

1976 

JUN 

19 

1025 

29 

7. 

1.71 

EXPLOSION 

1976 

AUG 

22 

1014 

30 

6. 

2.26 

EARTHQUAKE, 

AZ 

45,  MAG  4.0,  D  340 

1976 

AUG 

26 

1430 

31 

140. 

1.47 

EXPLOSION 

1976 

OEC 

08 

1460 

32 

69. 

1.12 

EXPLOSION 

1976 

DEC 

21 

1509 

33 

31. 

0.98 

EXPLOSION 

1977 

JAN 

13 

0709 

34 

13. 

2.00* 

EARTHQUAKE, 

AZ 

170,  MAG  3.8,  D  370 

1977 

JAN 

13 

2009 

35 

6. 

3.96 

EARTHQUAKE, 

AZ 

330,  MAG  3.7,  D  377 

1977 

FEB 

16 

1763 

36 

25. 

1.25 

EXPLOSION 

1977 

MAY 

31 

1640 

37 

3. 

3.67 

EARTHQUAKE, 

AZ 

330,  MAG  3.7,  D  330 

1977 

AUG 

04 

1640 

38 

130. 

0.86* 

EXPLOSION 

1977 

AUG 

12 

0220 

39 

40. 

3.73 

EARTHQUAKE, 

AZ 

170,  MAG  4.8,  D  421 

1977 

AUG 

19 

1765 

40 

471. 

1.27* 

EXPLOSION 

1977 

SEP 

15 

1436 

41 

18. 

1.28 

EXPLOSION 

1977 

OCT 

26 

1415 

42 

31. 

1.54 

EXPLOSION 

1977 

NOV 

10 

0235 

43 

14. 

3.28 

EARTHQUAKE, 

AZ 

126,  MAG  4.0,  D  302 

1977 

NOV 

17 

1930 

44 

40. 

1.18 

EXPLOSION 

1977 

DEC 

14 

1530 

45 

496. 

1.13 

EXPLOSION 

1978 

FEB 

14 

0435 

46 

16. 

2.96 

EARTHQUAKE , 

AZ 

45,  MAG  4.8,  0  340 

1978 

MAY 

23 

0647 

47 

20. 

4.45 

EARTHQUAKE, 

AZ 

20,  MAG  4.6,  D  434 

1978 

JUN 

16 

0421 

48 

64. 

2.17 

EARTHQUAKE, 

AZ 

170,  MAG  4.3.  D  345 

1978 

JUL 

07 

1400 

49 

24. 

1.16* 

EXPLOSION 

1978 

JUL 

29 

2232 

50 

52. 

2.08 

EARTHQUAKE , 

AZ 

45,  MAG  4.0,  D  340 

1978 

AUG 

01 

0902 

51 

40. 

3.87 

EARTHQUAKE, 

AZ 

350,  MAG  4.6,  D  397 

1978 

AUG 

01 

0947 

52 

24. 

5.27 

EARTHQUAKE , 

AZ 

350,  MAG  4.5,  D  397 

1978 

AUG 

01 

1026 

53 

16. 

2.21 

EARTHQUAKE , 

AZ 

350,  MAG  4.2,  D  397 

1978 

AUG 

31 

1400 

54 

1500. 

0.67 

EXPLOSION 

1978 

AUG 

31 

2356 

55 

21. 

0.50 

COLLAPSE 

1978 

NOV 

02 

1525 

56 

32. 

0.71 

EXPLOSION 

1979 

JAN 

06 

0120 

57 

16. 

2.55 

EARTHQUAKE, 

AZ 

75,  MAG  4.2,  D  377 

1979 

JAN 

24 

1800 

58 

25. 

1.39 

EXPLOSION 

1979 

FEB 

22 

0716 

59 

8. 

4.16 

EARTHQUAKE, 

AZ 

0,  MAG  3.5,  D  230 

1979 

FEB 

22 

1657 

60 

700. 

2.54 

EARTHQUAKE . 

AZ 

0,  MAG  5.0,  D  230 

*  -BY  CONVERSION  TO  JAMESTOWN  SHORT  PERIOD  (SEE  TEXT) 


Table  2.  Same  format  as  Table  1  for  the  S/P  discrimination. 


Figure  3.4.  This  figure  used  to  estimate  the  significance  of 
the  data  presented  in  Figure  3.3.  Where  three  readings  of  the 
ratio  of  S  to  P  amplitude  were  available,  a  standard  deviation 

was  computed  and  plotted  here.  The  horizontal  lines  delimit 
the  decision  band:  earthquakes  are  above  it  and  explosions  below. 
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measured,  while  In  Figure  3.1  we  are  required  to  calibrate  the  plot  for  each  re¬ 
gion  (which  Is  Impossible  In  the  regions  of  Interest);  (5)  the  phases  measured, 
l.e.,  Pg  and  Sg  (or  "Lg"  In  the  parlance  of  Russian  and  Midwestern  American  se¬ 
ismologists)  will  have  envelope  amplitudes  not  strongly  dependent  on  source 
focal  mechanism  as  they  are  crustal  channel  waves  that  sample  the  focal  sphere 
at  many  takeoff  angles.  The  simple  explanation  of  why  the  discriminant  works  Is 
In  the  energy  leaving  the  source  (hence  available  to  travel  outward  as  Pg  or 
Sg)«  shear  sources  send  out  far  more  (factor  of  ten)  energy  In  S  than  In  P. 
This  Is  the  same  reason  given  far  the  efficiency  of  the  Msimb  discriminant  as 
noted  above. 

It  Is  worth  noting  that  the  data  discussed  In  this  section  raises  serious 
questions  as  to  whether  we  understand  In  even  a  fundamental  way  the  nature  of 
surface-wave  excitation  by  explosions.  Although  It  Is  claimed  by  some  that  the 
S  leaving  the  source  excites  the  bulk  of  the  surface  waves  leaving  a  seismic 
source  (hence  the  Mssmb  discriminant),  I  saw  many  explosions  with  essentially  no 
Sg  phase  that  generated  huge  surface  waves  (hence  failure  of  the  discriminant  In 
Figure  3.1).  That  Is,  the  Sg  phase  seen  In  explosions  appears  uncorrelated  with 
the  amplitude  of  the  surface  waves.  I  believe  that  the  mechanism  of  surface 
wave  excitation  from  explosions  departs  slgnlf leant  ly  from  processes  which  can 
be  described  by  first-order  or  linear  theory  (which  encompasses  almost  the  to¬ 
tality  of  the  work  done  at  near-regional  distances  In  seismology).  As  I  have 
looked  at  more  of  this  data,  I  become  more  and  more  convinced  that  Vlecelll's 
(1973)  treatment  of  the  generation  of  these  waves  through  the  mechanism  of  spall 
closure  Is  a  better  route  toward  understanding  this  phenomenon.  A  similar  argu¬ 
ment  may  apply  also  to  the  generation  of  body  waves,  but  the  modelling  done  so 
far  seems  more  self-consistent  to  me  than  the  work  on  surface  waves.  Because  of 
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these  feelings,  I  would  recommend  stronger  emphasis  on  discrimination  methods 
which  utilize  body  waves,  as  I  believe  that  such  methods  are  more  likely  to  be 
capable  of  realistic  theoretical  treatment. 

Finally,  I  would  like  to  qualify  what  may  appear  to  be  an  enthusiastic  en¬ 
dorsement  of  the  S/P  method  of  dl scr Imlnatlon  discussed  above.  First,  the  popu¬ 
lations  do  not  separate  by  any  very  large  margin,  so  failures  of  the  discrimi¬ 
nant  are  bound  to  occur  In  other  areas  where  the  populations  move  slightly 
closer  together.  Second,  the  failure  to  discriminate  Event  28  In  Table  2  Is 
particularly  disturbing,  as  It  was  the  natural  earthquake  closest  of  any  studied 
here  to  the  test  site.  Bakun  and  Johnson  (1970)  had  found  a  dandy  discriminant 
that  turned  out  to  fall  for  events  on  the  test  site  (l.e.,  the  Massachusetts 
Mountain  earthquakes  were  Identified  as  explosions  by  their  method*  unpublished 
data,  1973)  as  were  those  aftershocks  of  nuclear  explosions  that  were  not  col¬ 
lapses.  Of  note  Is  that  Murphy  and  Lahoud  (1975),  using  a  rather  more  sophisti¬ 
cated  procedure  based  on  the  ratio  of  S  to  P,  clearly  dlscr Imlnated  the  Masse- 
chusetts  Mountain  malnshock  at  regional  distances  (their  plots  show  an  S  to  P 
ratio  of  4  or  more).  It  Is  probable  that  refinements  along  the  lines  presented 
by  these  authors  would  lead  to  a  s Ignlf leant ly  more  robust  discriminant. 

As  a  final  note  of  pessimism,  note  In  Figure  3.3  that  the  S/P  ratio  seems 
to  be  converging  at  the  small  magnitude  end  of  the  plot,  precisely  the  area  of 
greatest  interest. 

4.  SEISMIC  DISCRIMINATION  AT  REGIONAL  DISTANCES  USING  PN 


In  this  section  I  describe  work  aimed  at  obtaining  a  seismic  discriminant 
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based  on  analysis  of  the  PN  phase  alone.  The  study  Is  motivated  by  some  theo¬ 
retical  and  practical  considerations.  First,  PN  leaves  the  source  at  rather 
steep  angles  of  Incidence;  thus  we  sample  a  rather  confined  area  of  the  focal 
sphere,  which  has  the  potential  for  Increasing  the  precision  of  the  determina¬ 
tion  relative  tc  measurement  of,  say,  the  Pg  phase.  Second,  PN  travels  deep  In 
the  crust,  thus  (perhaps)  avoiding  some  of  the  uppei — crust  complexities  that  are 
sure  to  hinder  detailed  understandl ng  of  the  Pg  phase.  For  the  purpose  of  rou¬ 
tine  Identification  and  discrimination,  the  PN  phase  Is  an  attractive  one  to 
consider  and,  so  far  as  I  know,  such  a  study  has  not  previously  been  attempted. 

The  method  to  be  used  Is  Inversion  for  the  seismic  moment  tsnsor  as  des¬ 
cribed  in  previous  sections.  To  this  end  we  have  selected  two  low-  yield  events 
on  the  test  site,  one  <12  May  1976  at  1960  GCT)  on  Yucca  Flat  and  the  other  (02 
Nov  19~8,  1626  GCT)  on  Pahute  Mesa.  Data  from  the  Lawrence  Livermore  Laboratory 
fast  tapes  was  digitized  at  260  samples  per  second  and  decimated  to  60 
samples/sec  after  digital  alias  filtering  at  10  Hz.  Three-component  data  at 
Mina,  Kanab,  Landers,  and  Elko,  distributed  around  the  test  site  In  azimuths  at 
distances  of  200  to  400  km,  was  obtained.  This  data  provides  an  Ideal  proving 
ground  for  discriminants  of  this  kind,  as  the  data  quality  Is  excellent. 

The  discriminant  I  Intend  to  test  Is  based  on  the  prediction  of  ground  mo¬ 
tion  from  the  canonical  sources  required  for  Input  to  Brian  Stump's  Inversion 
cede  (Stump  and  Johnson,  1977).  Accordingly,  Green’s  functions  were  prepared 
using  MEXEC  for  a  model  similar  to  that  presented  by  Priestley  and  Brune  (1978) 
based  on  Great  Basin  surface  wave  analysis,  l.e.,  a  3-layer  crust  over  a  mantle 
with  a  PN  velocity  of  7.8  km/sec.  Due  to  anisotropy  or  some  other  failure  of 
the  Priestley  and  Brune  model  at  a  precision  of  tenths  of  a  second,  predicted 
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onset  times  of  PN  were  Inconsistent  at  the  four  stations.  Therefore,  I  adjusted 
the  times  so  that  the  onset  times  of  the  theoretical  and  observed  data  agreed  to 
avoid  spurious  offsets  In  phase  that  would  be  returned  by  moment  tensor  Inver* 
s  i  on . 


It  Is  clear  that  the  present  exercise  Is  going  to  be  a  severe  test  of  the 
moment  tensor  method.  Consider,  for  example,  Figures  4.1  and  4.2.  In  this  and 
other  figures,  groups  of  four  traces  are  shown.  The  top  trace  Is  theoretical 
ground  motion  at  the  appropriate  receiver  site  for  a  step  function  source.  The 
next  trace  Is  the  result  of  passing  the  theoretical  trace  trace  through  the  In¬ 
strument  (velocity  flat  systems  In  this  case).  The  third  trace  shows  the  effect 
of  Inserting  a  finite  source  time,  and  should  be  compared  with  the  fourth  trace, 
the  observations.  In  Figures  4.1  and  4.2  we  compare  vertical  components  for  the 
02  November  explosion  as  recorded  at  Mina  and  Kanab  (left  and  right  Images  of 
Figure  4.1)  and  Landers  and  Elko  (left  and  right  Images  of  Figure  4.2).  Note  In 
Figure  4.1  that  the  theoretical  trace  shown  Includes  a  part  of  Pg,  so  that  Pn  Is 
not  clearly  seen.  The  steppy  synthetic  ground  motions  on  the  other  three  sta¬ 
tions  are  what  Is  expected  by  the  theory  (PN  at  distance  mirrors  motion  at  the 
source ) .  Note  that  the  theoretical  PN  waveforms  are  all  similar,  but  the  data 
at  the  four  stations  are  quite  different.  Note  also  the  significant  difference 
In  frequency  content.  Theoretical  source  rise  times  of  0.4  second  and  overshoot 
ratios  of  about  1.5  to  1  were  Included  In  the  convolution  for  source  finiteness 
(third  traces  of  each  Image).  Numbers  given  to  the  right  of  each  trace  Indicate 
ma<1mum  trace  amplitudes  (the  bottom  trace  Is  scaled  to  ground  velocity  In  cgs 
units)!  note  the  variations  In  these  from  station  to  station  (bottom  traces). 
For  generalized  ray  theory  to  fit  these  observations,  a  crustal  model  with  con¬ 
siderably  finer  detail  will  have  to  be  employed. 
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l\SHrP  4^—  ^ame  format  as  Figure  4.1,  but  for  the  vertical  components  at  Landers  (left)  and 
Elko  (right).  Numbers  right  of  each  trace  are  maximum  trace  amplitude.  The  data  traces  are 
in  units  of  velocity  (pm/sec).  I 


F-l9Mre  4.3.  This  and  remaining  figures  of  this  section  are  in  the  same  format,  vertical 
components  on  the  left  and  radial  components  on  the  right.  Traces  from  top  to  bottom  are 
as  Figures  4.1  and  4.2,  i.e.,  theoretical  traces  on  top  and  the  data  on  the  bottom. 

Shown  for  the  Landers  station  are  the  data  and  the  theoretical  displacements  to  be  expected 
from  a  strikeslip  earthquake  oriented  45  degrees  off  the  line  from  station  to  epicenter. 


Figure  4.4.  Same  as  Figure  4.3,  but  for  the  radial  and  vertical  components  of  a  dipslip  earthquake 
with  fault  plane  oriented  perpendicular  to  the  line  from  the  station  to  epicenter.  Landers 
station. 
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Figure  4.5.  Same  format  as  Figure  4.3,  again  at  Landers,  but  for  a  source  consisting  of  a 
compensated  linear  vector  dipole,  i.e.,  gil.l  +  gi2,2  -  2*gi3,3. 
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In  the  set  of  figures  4.3  to  4.6.  we  compare  the  appearance  of  the  PN  phase 
produced  by  explosion  and  earthquake  sources  at  Landers  for  the  02  November  ex¬ 
plosion.  In  each  figure  we  present  the  vertical  component  on  the  left  and  the 
radial  component  on  the  right.  Figure  4.2  Is  for  a  strlkesllp  earthquake  ori¬ 
ented  at  46  degrees  from  the  receiver.  Figure  4.3  Is  for  a  dlpsllp  earthquake 
whose  fault  plane  Is  oriented  In  line  with  the  receiver.  Figure  4.4  Is  the  mo¬ 
tion  from  a  compensated  linear  vector  dipole.  Finally,  Figure  4.6  Is  for  the 
explosion  source.  Note  that  the  theoretical  records  for  earthquake  sources  are 
considerably  more  complex  than  those  for  the  explosion,  mainly  because  the 
latter  Includes  no  generalized  raypaths  that  leave  the  source  as  an  SV  wave. 
However,  on  all  these  sources  the  motions  In  the  PN  phase  observed  lasts  consid¬ 
erably  longer  than  the  theoretical  calculation  Indicates,  and  none  of  them  looks 
particularly  like  the  data. 

All  of  the  data  and  Green's  functions  shown  In  this  section  have  been 
placed  on  the  Cyber  176  computer  at  Klrtland  Air  Force  Base  In  Albuquerque,  and 
await  the  running  of  Brian  Stump's  code.  In  a  subsequent  report  the  outcome  of 
this  and  any  follow-on  experiments  will  be  described. 

S.  ANALYSIS  OF  CLOSE-IN  ACCELEROMETER  RECORDS  OF  MEGATON  EXPLOSIONS 

Peppln  (1977),  Stump  (1979),  and  Helmberger  and  Hadley  (1900)  have  present¬ 
ed  analyses  of  close-in  (8  km)  records  of  the  megaton  explosions  JORUM  and  HAND- 
LEY,  which  occurred,  respectively  on  16  Sep  1969  at  1430  GCT  and  23  Mar  1970  at 
1900  GCT.  These  analyses  have  grown  ever  more  sophisticated.  Peppln's  solution 
Is  the  simplest.  That  of  Stump  provides  the  best  fit  of  theory  with  data; 
Helmberger  and  Hadley  have  provided  the  most  realistic  Green's  functions. 
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Results  of  these  studies  are  In  significant  dl sagreement .  While  the  latter 
study  concludes  that  overshoot  at  the  source  Is  necessary  to  fit  the  observa¬ 
tions,  the  former  two  do  not.  The  problem  Is  well  worth  Investigating,  because 
tills  overshoot  Is  relevant  to  the  under  stand  1 ng  of  body-wave  excitation  by  ex¬ 
plosions.  In  this  section  I  describe  efforts  aimed  at  merging  the  calculations 
of  realistic  Green's  functions  using  the  MEXEC  code  with  the  powerful  moment 
tensor  Inversion  method  Implemented  by  Stump.  Material  In  this  section  Is  pre¬ 
sented  at  the  American  Geophysical  Union  meetings  In  San  Francisco  (Peppln, 
1979). 

In  Section  1  I  described  the  development  of  the  MEXEC  code.  This  code  per¬ 
mits  the  determination  of  multilayer  Green's  functions  at  any  distance  range 
(although  this  Is  not  possible  as  a  practical  matter  for  reasons  of  numerical 
stability  and  cost  In  a  wide  variety  of  Interesting  problems).  It  was  specifi¬ 
cally  written  for  the  problem  described  in  this  section.  We  have  computed  four 
test  cases *  source  In  a  layer,  source  below  the  layer,  Helmberger-Hadley  model, 
and  extended  He lmberger-Had ley  model.  These  are  discussed  In  subsequent  subsec¬ 
tions. 


5.1.  Source  In  a  Layer 

Iri  Figure  5.1  we  summarize  the  layer  models  used  In  the  various  test  cases. 
For  this  section,  we  take  the  model  In  the  left-center  part  of  the  figure,  with 
the  original  halfspace  model  used  by  Peppln  <1977)  above  for  comparison.  Model 
parameters  are  arbitrarily  chosen  so  that  the  onset  time  Is  as  observed  on  the 
test  site,  and  with  Poisson's  ratio  1/4,  In  Figure  5.2  we  show  the  comparison 
of  the  data  with  the  Green's  function  computed,  which  contains  all  first-order 
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Figure  5.1 .  Layer  models  considered  in  this  section  drawn  with  layer  thicknesses  to  scale. 
Stars,  represent  source.  Top  model:  halfspace  case  considered  by  Peppin  (]977);  other  moc 
as  discussed  in  the  text.  Triangles  show  locations  of  extended  Helmberger-Hadley  model. 
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reflections  In  the  layer  Including  the  free-surface  reflection.  As  In  the  pre¬ 
vious  section,  vertical  data  Is  shown  on  the  left,  radial  on  the  right.  The  top 
trace  Is  theoretical  response  to  a  step  function;  the  second  trace  results  from 
passing  the  first  trace  through  the  accelerometer 5  the  third  trace  shows  the 
effects  of  a  finite  source  duration;  the  fourth  trace  Is  data  taken  on  the  test 
site. 


Note  that  a  rise  time  of  0.6  second  has  been  used  In  Figure  5.2,  and  an 
overshoot  ratio  of  1.06  to  1  (almost  no  overshoot).  If  more  overshoot  Is  added 
the  fit  to  the  first  second  of  the  data  Is  degraded  (Figure  5.3).  Although  the 
fit  on  the  waveforms  of  the  P  onset  Is  pretty  good,  the  relative  ratios  predict¬ 
ed  fall  to  match  the  observed  by  a  large  factor  (note  how  the  theoretical  radial 
Is  larger,  but  the  observed  vertical  data  Is  larger*  numbers  right  of  etch 
trace  in  these  figures).  Thus,  although  overshoot  at  the  source  Is  evidently 
not  necessary  here,  the  ratio  of  vertical  to  radial  motion  provides  a  serious 
problem,  the  same  one  discussed  by  Peppln  (1977). 

5.2.  Source  Below  a  Single  Layer 

In  Figure  5.4  we  have  placed  the  source  below  a  slow  layer*  see  the  layer 
model  center  right  In  Figure  5.1.  Layer  parameters  are  chosen  so  that*  (1)  the 
P  onset  times  agree  with  the  observed,  (2)  the  vertical  component  of  P  In  the 
synthetics  has  about  the  correct  amplitude  relative  to  the  radial,  and  (3)  Pois¬ 
son's  ratio  Is  1/4.  The  Green's  functions  Include  all  directs  and  first-order 
multiples  (with  conversions)  In  the  top  layer.  Note  the  radial  component  In 
Figure  5.4.  The  theoretical  trace  fits  the  data  as  well  as  any  In  this  section 
right  through  the  Rayleigh  wave  arrival  In  character,  and  In  fair  detail  near 
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Figure  5.3.  Same  format  as  Figure  5.2  and  same  layer  structure,  but  now  the  source  has 
a  1.5  to  1  overshoot.  Note  that  the  P  onset  overshoots  too  much  compared  with  Figure 
5.2,  which  is  a  better  fit. 


Figure  5.1.  Same  format  as  Figure  5.2,  but  for  the  case  of  the  source  below  the  layer(right 
center  model  of  Figure  5.1).  Rise  time  of  0.6  second  as  before,  and  essentially  no  overshoot 
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the  P  onset.  However,  the  vertical  component  synthetic  looks  only  remotely  like 
the  data.  The  difference  in  signal  character  between  the  two  synthetics  is 
(probably  artificially)  caused  by  the  sense  of  the  P  to  S  conversion;  for  the 
radial,  this  comes  in  about  right,  but  for  the  vertical  it  is  Inverted  relative 
to  P  (second  spike  of  the  top  trace  in  Figure  6.4).  Adding  overshoot  at  the 
source  did  not  help  to  improve  the  fit,  nor  did  varying  the  rise  time.  What  is 
shown  Is  the  best  I  could  find  for  this  layer  model.  Note  the  we  11 -pronounced 
near-field  effects  on  the  top  traces  (the  long-period  trends  in  the  trace  before 
the  surface  wave). 

5.3.  Helmberger-Hadley  model 

I  computed  Green's  functions  for  the  Helmberger-Hadley  (1980)  model  of  the 
test  site  (Figure  6.1  bottom),  summing  directs,  first-order  reflections,  and 
free-surface  reflections.  The  resulting  Green's  functions  can  be  seen  in  Figure 
6.6.  As  can  be  seen  the  fit  Is  quite  poor,  and  probably  because  our  source  Is 
placed  in  the  third  layer  rather  than  the  second.  We  would  have  hoped  for  a 
better  comparison  with  Helmberger  and  Hadley,  and  the  failure  to  attain  this  may 
be  causod  byi  (1)  the  fact  that  our  fit  is  done  In  acceleration  and  theirs  is 
in  velocity,  (2)  the  different  shape  of  our  source  time  function,  or  (3)  an 
error  In  the  computation  of  the  Green's  function.  During  the  next  contract  per¬ 
iod  I  will  conduct  some  definitive  checks  with  existing  genera  1 1  zed-ray  codes 
from  other  groups,  although  I  doubt  that  the  MEXEC  code  can  be  grossly  in  error 
at  this  point. 


5.4.  Extended  Helmberger-Hadley  model 


H|ure  5^f.  Same  format  as  Figure  5.2  for  the  Helmberger-Hadley  model  (lower  one  in  Figure 
5.1 The  source  has  mild  (1.05  to  l)overshoot  and  a  rise  time  of  0.6  second. 
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It  Is  known  (e.g.,  Rodean,  1971)  that  tha  sourca  of  a  nuclear  explosion  Is 
no  way  a  point  source.  Thus,  all  of  the  above  studies  which  have  employed  this 
In  tha  modal  violate  known  physics.  In  tha  previous  subsections  we  have  seen 
how  critical  tha  appearance  of  the  Green's  functions  depends  on  the  placement  of 
the  source,  whether  In  the  top,  second,  or  third  layer.  To  Insure  that  the 
Green's  functions  are  not  dominated  by  this  (arbitrary)  placement  In  the  layer 
stack,  we  must  find  some  rational  way  to  account  for  the  source  finiteness. 

I  have  here  used  a  crude  but  reasonable  approach.  Suppose  a  source  radius 
of  0.6  km  following  Peppln  (1977).  Then  the  source  will  extend  from  the  top  of 
the  second  layer  to  the  bottom  part  of  the  third  one  (see  bottom  Image  of  Figure 
6.1).  Thus,  replace  the  point  source  by  a  pair  of  point  sources,  one  midway 
through  the  second  layer  and  one  at  a  depth  of  1.63  km  In  the  third  layer. 
Scale  these  sources  by  the  product  of  the  volume  In  the  respective  layers  and 
the  rigidity  of  the  layers,  so  that  27%  Is  In  the  second  layer  and  73%  Is  In  the 
bottom  layer  (deltas  In  Figure  6.1  show  the  locations  of  these  sources). 
Results  are  presented  In  Figure  6.7.  The  fit  Is  still  poor,  but  a  little  better 
than  the  fit  obtained  In  the  previous  section. 

It  seems  clear  that  the  synthetics  are  most  like  the  data  when  the  source 
Is  placed  In  the  second  (not  top  or  third)  layer.  In  contrast  to  the  statements 
of  Peppln  (1977),  It  Is  also  clear  that  we  need  not  resort  to  compound  sources 
to  fit  these  observations,  an  explosion  will  do  very  nicely.  Of  note  Is  that 
the  synthetics  for  earthquakes  look  nothing  like  the  data  for  any  of  these  mo¬ 
dels;  Figure  5.8  Is  typical.  Here  we  show  the  response  to  be  expected  from  a 
dlpsllp  earthquake  placed  at  the  same  source  depth  as  the  explosions.  Energy  In 
the  S  and  surface  waves  dominates  the  synthetics,  and  In  particular  dominates 
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the  P  arrival.  It  Is  certain  that  moment  tensor  Inversion  on  any  of  these 
Green's  functions  will  result  In  a  source  dominated  by  the  Isotropic  (l.e.,  ex¬ 
plosive)  part.  That  Is  the  most  Important  question  to  explore  In  evaluating  the 
moment  tensor  Inversion  method. 

All  of  the  Green's  functions  and  data  shown  In  this  section  are  on  the 
Cyber  176  computer  at  Klrtland  Air  Force  Base  In  Albuquerque  awaiting  Input  to 
Brian  Stump's  moment  tensor  Inversion  code. 

6.  PN  DISCRIMINATION  IN  A  SHIELDLIKE  ENVIRONMENT 

An  effort  has  been  made  during  this  contract  period,  through  phone  communi¬ 
cations  and  a  visit,  to  obtain  data  from  St  Louis  University  on  the  deep  Merrl- 
mac  mine  blasts  of  northern  Missouri.  These  provide  some  of  the  most  provoca¬ 
tive  data  for  seismic  discrimination  In  existence,  as  the  blasts  are  large 
(fraction  of  a  klloton)  and  recorded  clearly  out  to  ranges  of  hundreds  of  km. 
Records  of  these  blasts,  fired  3  km  underground,  are  striking  In  the  large 
S-waves  present;  they  would  be  classified  as  explosions  using  the  S/P  discrimi¬ 
nant  discussed  In  Section  3,  looking  very  much  like  local  earthquakes  Indeed. 
Thus,  this  data  should  provide  another  severe  test  of  the  moment  tensor  analysis 
method.  PN  data  Is  available  on  a  dozen  or  so  stations  to  the  south  and  east  of 
these  blasts,  as  well  as  for  some  nearby  earthquakes  f»r  comparison.  I  have  at¬ 
tempted  to  obtain  this  data  for  Investigation  of  PN  discrimination  as  described 
In  Section  4  of  this  report. 

Due  to  limitations  In  staff  and  equipment  at  St  Louis,  the  logistics  of 
this  effort  are  complicated  and  fairly  expensive.  We  must  rent  dual  tape  play- 
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back  systems  at  considerable  cost  and  have  these  shipped  to  St  Louis  for 
transfer  of  the  analog  data  to  dub  tapes.  This  will  add  a  significant  expense 
to  this  contract,  of  the  order  $2,000.  The  data  will  be  digitized  and  processed 
here  on  the  Selsmologlcal  Laboratory  computer  system  sometime  after  the  first  of 
the  year. 

7.  PERIPHERAL  STUOIES  OF  INTEREST  TO  THIS  CONTRACT 

In  this  section  Is  described  work  that  Is  off  the  major  emphasis  of  this 
contract,  but  nevertheless  Is  relevant.  We  describe  efforts  aimed  at  acquisi¬ 
tion  of  high-quality  close-in  data  for  earthquakes,  a  problem  on  earthquake 
spectral  corner  frequencies,  and  the  use  of  explosions  for  the  determination  of 
phase  velocity  on  the  test  site. 

7.1.  Digital  data  Acquisition  and  Spectral  Corner  Frequencies 

Over  the  p_  t  year  and  a  half,  1  have,  when  time  permits,  attempted  to  con¬ 
tinue  our  program  of  data  acquisition  using  the  high-quality  digital  event  re¬ 
corders  built  in  1976  and  1977  with  support  from  AFOSR.  These  Instruments  have 
by  now  gone  far  In  meeting  the  objectives  of  the  Near-Field  Project  sponsored  by 
the  Air  Force  In  the  early  70s.  We  now  have  about  1,000  records  of  earthquakes 
In  the  Sierra  and  Great  Basin  region,  many  of  unprecedented  quality.  This  ef¬ 
fort  has  become  timely  In  view  of  the  significant  sequence  of  earthquakes  which 
have  been  occurring  In  the  Sierra  Nevada  boundary  zone.  The  records  obtained 
Include  many  when  the  system  was  running  flat  In  seismic  displacement  from  10 
seconds  to  50  Hz,  extremely  wide  band  for  a  portable  system.  The  system  has 
proven  capable  of  fulfilling  the  design  requirements  of  the  Near-Field  Project; 
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an  earthquake  of  magnitude  3.3  recorded  directly  under  the  recorder  was  within 
full  scale  by  a  factor  of  6j  see  Figure  7.1. 

The  acquisition  has  lead  to  a  number  of  papers  including  one  involving  this 
author  (Somerville  et  al.,  1980).  In  this  paper  we  discuss  selsmotectonics  of 
the  conspicuous  Genoa-Carson  fault  system.  In  response  to  the  occurrence  of  a 
sequence  on  the  south  end  of  this  zone,  we  placed  a  digital  event  recorder  al¬ 
most  directly  over  the  source  of  the  earthquakes,  which  was  found  to  be  a  tight 
cluster  2  km  In  extent  and  10  km  deep.  In  Figure  7.2  we  show  a  typical  record 
obtained  for  an  earthquake  of  about  magnitude  1.5.  Note  the  extreme  complexity 
of  the  records  for  the  nearly  vertical  wave  propagation. 

The  major  relevance  of  this  work  to  source  theory  (hence,  e.g.,  better 
understanding  of  seismic  discrimination),  Is  In  further  documentation  of  the 
rather  enigmatic  behavior  of  the  ratio  *R“  of  P  to  S-wave  spectral  corner  fre¬ 
quencies.  Shown  In  Figure  7.3  are  the  results  of  the  Diamond  Valley  study, 
showing  large  <2. 6)  values  for  this  ratio.  Of  Interest  Is  that  this  seems  typi¬ 
cal  for  other  earthquake  sequences  recorded  In  the  Sierra  boundary  zone.  The 
result  Is  not  to  be  expected.  Peppln  and  Slmlla  (1976)  studied  earthquakes  from 
the  same  region  at  the  wideband  Jamestown  station,  which  entails  propagation 
across  the  Slerran  bathollth,  presumably  a  path  with  minimal  attenuation.  They 
found  nearly  unit  values  for  R.  That  Is,  more  P-wave  energy  Is  ev’dent  at  the 
close-in  stations  than  Is  evident  at  the  more  distant  stations. 

It  seems  unlikely  that  this  can  be  the  result  of  anelastlc  attenuation, 
since  most  work  Indicates  that  S  energy  attenuates  as  least  as  rapidly  as  P  en¬ 
ergy.  If  It  Is  a  source  effect,  then  It  seems  puzzling  that  we  again  see  about 
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unit  values  for  the  ratio  R  In  the  case  events  recorded  at  comparable  distances 
near  The  Geysers,  California*  see  Figure  7.4  from  Peppln  and  Bufe  (1980),  using 
Identical  Instrumentation  and  processing  as  was  used  to  obtain  Figure  7.2  for 
the  Diamond  Valley  earthquakes. 

In  recent  source  theories  the  ratio  R  proves  to  be  a  strong  diagnostic  as 
to  the  nature  of  rupture  propagation  at  the  source.  The  results  we  have  pre¬ 
sented  here,  and  others  we  are  preparing  for  five  other  sequences  In  the  Sierra 
province,  may  be  leading  to  an  Important  result  either  In  terms  of  knowledge  of 
the  source  or  In  terms  of  wave  propagation  and  attenuation  In  the  lower  crust 
and  upper  mantle  of  the  Basin  and  Range. 

7.2.  Analysis  of  Explosion  Surface  Waves  for  Phase  Velocity 

The  existence  of  ultralong  period  data  at  the  Berleley  array  of  wideband 
seismometers  permits  Investigation  of  phase  velocity  over  short  paths  on  the 
test  site.  To  complement  the  work  of  Priestley  and  Brune  on  surface  wave  ana¬ 
lysis  In  the  Basin  and  Range,  I  have  studied  some  of  the  large  explosions  re¬ 
corded  at  Berkeley  and  Jamestown.  Data  quality  Is  excellent,  showing  clear  en¬ 
ergy  In  the  surface  wave  spectrum  of  explosions  out  to  periods  of  100  seconds. 
An  idea  of  the  data  quality  can  be  obtained  by  reference  to  Figure  7.6,  showing 
records  of  the  COLBY  test  taken  on  the  wideband  system  at  Jamestown.  Note  fac- 
tor-af-ten  signal  power  above  the  noise  out  to  frequencies  of  60  Hz. 

The  Idea  of  the  method  employed  here  was  to  use  two  explosions  on  the  test 
site  as  recorded  by  a  single  station  (Berkeley  or  Jame.town).  The  difference  In 
phase  between  similar  surface  wave  packets  of  explosions  along  a  line  of  azimuth 


Figure  7,j.  Ratio  of  spectral  corner  frequencies  of  P  and  S  waves 


Figure  7.<.  Vertical -component  surface  waves  for 
the  COLBY  test  as  recorded  on  the  Jamestown  displace 
"lent  transducer,  data  shown  above.  A  spectrum  with 
a  hundred-second  window  was  computed  below,  with  a 
noise  sample  shown  in  lighter  line  weight. 
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Figure  7.1.  Group  velocity  of  30  to  50-second  surface  waves  at  Jamestown  for 
the  HUtilSTER  test. 
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Figure  7.1.  Three  shots  on  almost  exactly  the  same  azimuth  line  from 
the  test  site  to  Jamestown,  with  trace  separation  proportional  to  the 
separation  of  the  respective  shotpoints.  Visual  correlations  used  to 
compute  phaee  velocity  in  the  4.2  to  8.6  second  range. 
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Figure  7.J.  Same  format  as  Figure  7.6  for  ' 
the  large  tests  MEENSTER  and  COLBY  as  recorded 
on  the  Berkeley  ultralong  period  system,  vertical 
component.  Periods  out  to  100  seconds  can  be 
seen  easily  in  the  spectral  records  of  these 
data. 
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gives  rise  to  a  measurement  of  the  structure  underneath  the  zone  between  the  two 
paths,  In  this  case  extending  to  the  bottom  of  the  crust  or  even  deeper.  While 
phase  resolution  Is  reduced  because  of  the  short  paths  used,  we  gain  the  advan¬ 
tage  of  precise  determination  of  phase  velocity  (In  theory)  for  areas  of  inter¬ 
est  on  the  test  site.  Group  velocity  analysis  was  attempted  first  on  the  re¬ 
cords  to  see  how  well  the  signals  were  dispersed!  see  a  typical  result  In  Fig¬ 
ure  7.6.  For  the  MUENSTER  test,  the  Information  extends  to  60  seconds. 

The  program  to  Invert  for  phase  velocity  was  not  successful.  However,  fn 
Figures  7.7  and  7.8  the  phase  moveout  can  clearly  be  seen  for  combinations  of 
explosions;  the  surface  wave  packets  of  the  different  explosions  look  very  sim¬ 
ilar  Indeed,  and  phase  correlation  Is  easy  to  do  visually.  Our  aim  Is  to  extend 
and  overlap  the  extensive  data  set  of  Priestley  and  Brune,  hopefully  to  40  sec¬ 
onds,  with  emphasis  for  paths  on  the  test  site.  A  cooperative  effort  with  a 
graduate  student  at  Berkeley  was  1nT<ated  during  this  contract  period  for  the 
purpose  of  pursuing  this  effort. 
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